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We consider colour-electric screening as expressed by the quark contribution to the Debye mass 
calculated in a PNJL model with emphasis on confining and chiral symmetry breaking effects. We 
observe that the screening mass is entirely determined by the nonperturbative quark distribution 
function and temperature dependent QCD running coupling. The role of the gluon background 
(Polyakov loop) is to provide strong suppression of the number of charge carriers below the tran¬ 
sition temperature, as an effect of confinement, while the temperature dependent dynamical quark 
mass contributes additional suppression, as an effect of chiral symmetry breaking. An alternative 
derivation of this result from a modified kinetic theory is given, which allows for a slight generaliza¬ 
tion and explicit contact with perturbative QCD. This gives the possibility to gain insights into the 
colour screening mechanism in the region near the QCD pseudocritical temperature and to provide 
a guideline for the interpretation of lattice QCD data. 

PACS numbers: 12.38.Mh,12.38.Gc 


I. INTRODUCTION 


Theoretical and experimental investigations of quan¬ 
tum chromodynamics (QCD) at finite temperatures are 
performed with the aim to gain insights to the mecha¬ 
nisms of chiral symmetry restoration and deconfinement. 
From this perspective the heavy-quark (HQ) potential 
is one of the most important probes to be studied. In 
the vacuum the heavy quark-antiquark system is well de¬ 
scribed in terms of effective field theories, owing to the 
energy scale separation mm and the fact that the HQ 
potential can be properly defined in terms of large Wil¬ 
son loops. Then the spectrum of charmonium and bot- 
tomonium states as solutions of the Schrodinger equation 
for heavy quarkonia can be extracted and faced experi¬ 
ment. Finite temperature studies rely on the extraction 
of a static potential from ab initio simulations of lat¬ 
tice QCD (IQCD) considering the singlet free energy due 
to a pair of static color charges as a function of their 
distance by means of Polyakov loop correlators [S] . How¬ 
ever, the systematic field theoretic description within the 
NRQCD framework [2] is much more subtle and delicate 
because new variable energy scales like the temperature 
T or the screening mass enter the problem distorting 
the energy hierarchy. In addition, from the definition of 
a real-time potential [1] it was realized that also an imag¬ 
inary part in the potential appears at finite temperature 
[SHU. Using lattice QCD studies of this complex values 
static potential it was recently found that the real part is 
well described by the color singlet free energy of a static 
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quark-antiquark pair [5]. 

While those results shed light on the important ques¬ 
tion about the identification of the HQ potential with the 
colour singlet free energy another important question 
is that of microphysics insights into the screening mecha¬ 
nism. Based on some first principle motivations |10j , Riek 
and Rapp have proposed an ansatz for the HQ potential 
m in the form of a screened Cornell potential where the 
Coulomb and the linear parts are subject to two differ¬ 
ent screening masses rriD and fhu, respectively. The fit 
of the temperature dependence of these parameters pro¬ 
vided in Ref. m using available IQCD data has revealed 
some unexpected aspects. First, the coulombic Debye 
mass mo has a linear behaviour with very small slope 
(smaller than expected from pQCD). Second, the screen¬ 
ing mass of the confining part, m_D, shows a strong sup¬ 
pression for temperatures below and a linear rise for 
high temperatures (higher than expected from pQCD). 

This is somehow different from the standard approach, 
where the IQCD data for the large distance part of 
this HQ potential are fitted either to a Debye screened 
Coulomb ansatz or to a form motivated by a Debye- 
Hiickel theory dSHIl]. An ansatz taking into account 
real and imaginary parts of the HQ potential has been 
recently considered in m in the context of quenched lat¬ 
tice data. Here we use moiT) from IQCD data with dy¬ 
namical fermions obtained from fitting the heavy quark- 
antiquark free energy at large separations to the standard 
Debye screened potential [T2] to compare with the model 
predictions. While its behaviour well above the pseu¬ 
docritical temperature of the QCD phase transition 
can qualitatively be understood in terms of perturbation 
theory, the interpretation of the lattice data in the vicin¬ 
ity of Tc require essentially nonperturbative approaches 
addressing effects of confinement and chiral symmetry 
breaking. The leading order perturbative result reads 
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mo = (1 + while the next-to-leading order can 

be obtained by resummation of the leading contribution 
of the high temperature expansion [Mil]- A detailed 
understanding of the physics behind the Debye mass in 
the nonperturbative domain is subject to many current 
studies. 

Approaches based on the operator product expansion 
(OPE) [IH], gauge/gravity duality [MED] or phenomeno¬ 
logical models uni make an attempt to give a microscopic 
description of the screening phenomenon. However, the 
very definition and numerical determination of a screen¬ 
ing mass is obscured by the complications of the non- 
abelian nature of QCD and the large value of the coupling 
constant near the QCD transition region [211 Eg. 

In this paper we investigate screening effects in a PNJL 
model which proved successful in reproducing various as¬ 
pects of hadronic excitations in the medium |23| and of 
IQCD thermodynamics [24]. We will evaluate the one- 
loop polarization function using PNJL propagators with 
QED like vertices thus extending a previous calculation 
made for massless quarks [25] . In this rough way we im¬ 
plement confinement and chiral symmetry breaking ef¬ 
fects which in turn allows a comparison to the lattice 
data for the Debye mass. We will show that one can re¬ 
produce the correct shape of its temperature dependence. 
However, due to the absence of dynamical gluons in our 
PNJL model calculation, we lack dynamical degrees of 
freedom and therefore stay below the IQCD result. 

It is well known that the screening mass in QED or 
pQCD can also be derived from kinetic theory ESlEZl- 
Interrelations between plasma physics and quark gluon 
plasma are known to bring many relevant insights and 
physical motivation behind the held theory calculations 
[D5] . We also explore this possibility and modify the stan¬ 
dard kinetic theory approach by replacing usual Fermi- 
Dirac distribution functions for quarks with those mod- 
ihed by the coupling to the Polyakov loop. In this way 
we are able to reproduce the result for the Debye mass 
calculated within our model and furthermore to achieve 
contact with QCD by inclusion of effects of perturbative 
non-Abelian vertices. 

The paper is organized as follows. In section |T^ we 
outline the model calculation of the Debye mass within 
the PNJL model, whereby details are referred to the Ap¬ 
pendix. Section m presents the kinetic theory approach 


to the problem where we give a simple and intuitive 
derivation of the screening mass suggesting a straight¬ 
forward modification of the standard approach. Section 
|IV| is devoted to a comparison with IQCD data and their 
interpretation while section [Vj gives the conclusions. 


II. PNJL MODEL CALCULATION OF THE 
DEBYE MASS 

In reference [25] a model was considered where the vac¬ 
uum HQ potential was screened by a quark loop with 
internal lines coupled to a temporal background gluon 
field. For the static interaction potential V{q), = |qp, 

the statically screened potential is given by a resumma¬ 
tion of one-particle irreducible diagrams (RPA ’’bubble” 
resummation) 

Vsc{q) = V{q)/[l + F{0;q)/q^] , (1) 


where the longitudinal polarization function in the finite 
temperature case is defined via the projector decomposi¬ 
tion of the self energy (in Euclidean space) 


n^^(ia;, g) = F{iuj, 

+ G(ia;, 9 )Pj,, , 

( 2 ) 

where the projectors satisfy 


(pL )2 ^ ^ ^ 

pTpL ^ pLpT ^ Q ^ 

(3) 

+ P^v 

- 0 ^. ^2 ’ 

(4) 

jjT X QiQj 

Pj - Q2 

II 

II 

II 

0 

(5) 


Here, Q = (w, q) is the Euclidean four momentum. From 
this we get the gauge invariant longitudinal component 

F{iuj,q) = ^n 44 (ia;,q) , ( 6 ) 

and it can be calculated within thermal held theory as 


n44(*Wi;q) = ^ 


(Pp 

(2tt\ 


3 Tr[745'4,(ja;„;p)74S'$(iw„ - zw/ip - q)] . 


(7) 


Here w/ = 2'kIT are the bosonic and = (2n-|- l) 7 rT the 
fermionic Matsubara frequencies of the imaginary-time 
formalism. The symbol Tr stands for traces in color, ha- 
vor and Dirac spaces. 5'$ is the fermionic quark propaga¬ 


tor coupled to the homogeneous static gluon background 
held (^ 3 . Its inverse is given by [23JE2I 

S'J^(p;ia;„) = 7 • p - 740 ;^-I- 74 A 3 (/ 33 -I-m 


( 8 ) 
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where { 7 ^, 7 ;^} = —‘25^1, and ra = m{T) is the dynam¬ 
ically generated temperature dependent mass for light 
quarks as described, e.g., within the NJL model [3T]. The 
variable 1^3 is related to the Polyakov loop variable de¬ 
fined by [ 21 ] 

<i>(T) = = 1(1 + 2cos(/?^3)) , 


The physics of $(r) is governed by the temperature- 
dependent Polyakov loop potential U{^), which is fitted 
to describe the lattice data for the pressure of the pure 
glue system [21]. After performing the color-, flavor- and 
Dirac traces and making the fermionic Matsubara sum¬ 
mation we obtain [32] (see Appendix for the details) 


1144 (fw,g) 


rp^dp2MEp) 

9 / 2 z? 

Jo 7^ Ep 

j 4:Epiuj + + u)'^ — 4:Ep R_^_ (w) 1 

(9) 


where 



FIG. 1: Debye mass in leading order perturbation theory (dot¬ 
ted line), including the effect of dynamical chiral symmetry 
breaking without coupling to the Polyakov loop (dash-dotted 
line) and including the coupling to the Polyakov loop (solid 
line). Calculated with Us = 0.471 fitted to charmonium spec¬ 
trum at T = 0. 


i?±(a;) = — — 2iujEp ± 2pq , (10) 

and Re/(w) := i(/(a;)-I-/(—w)). Taking the static, long 
wavelength limit [23] [321133] we identify, after continua¬ 
tion fw —>■ go + to the Minkowski space, F(go = 0, g —>■ 
0 ) = m\){T) and whence 


The above result can be also obtained in a different 
way. It is well known that in QED the Debye mass is 
related to the pressure via the second derivative [32] [33] 

2 (rr \ 2^^E’(T, flf.) 

mD{T,Pe) = e -- , 13) 

dpi 


‘■D 


r dp [p^ + E^p] 
^ Jo 


MEp) 

E„ 


( 11 ) 


Here moiT) is the Debye mass, Ep = _|_ ^2 jg 

the quasiparticle dispersion relation for light quarks with 
Nc = 3 colour and Nf = 2 flavour degrees of freedom; 
f^{E) is the quark distribution function [33] 


ME) 


$(1 -b 2 e-^^)e-^^ -b e-^dE 
^1 -b 3$(1 -b e~P^)e~dE -p g-3/3E 


( 12 ) 


This result is different from the standard QED case only 
in that the Fermi-Dirac distribution has been replaced by 
the function (12). In comparison to the free fermion case 


PI33] the coupling to the Polyakov loop variable $(T) 
gives rise to a modification of the Debye mass, encoded in 
the modification of usual Fermi-Dirac distribution func¬ 
tion. In the limit of deconfinement ($ = 1), the case 
of a massive quark gas is obtained, while for confine¬ 
ment ($ = 0) one finds a considerable suppression. The 
temperature dependence of the Debye mass is shown in 
figure [^ and as expected from the very beginning it turns 
out to be below the free massless case, in the confined 
and in the transition region (with a pseudocritical tem¬ 
perature Tc « 200 MeV). For temperatures T » the 
free gas behaviour m'jy = 2/3g^T^ is reproduced. 


where here Pe is related to the electric charge of the 
system. This relation is a consequence of the Dyson- 
Schwinger equation for the photon self energy and the 
Ward identity relating the electron-photon vertex to the 
quark propagator. This is only true in abelian gauge the¬ 
ory and breaks down for non-abelian theories like QCD. 
On the technical level the proper Ward identity (called 
then Slavnov-Taylor identity) becomes much more com¬ 
plex and does not allow a simple derivation. Because 
our calculation in this section is similar to a one-loop 
QED calculation it is interesting to see if a similar re¬ 
lation holds also in the PNJL model. We check this in 
the finite temperature and zero chemical potential case 
by directly evaluating right hand side of eq. (13). Let us 
recall the quark contribution to the meanfield pressure 
reads |23] 


p,(r,rt = -2N,T j 

{In [1 -b 3($ -b $X_)A:_ -b xi] 

-b In [1 -b 3($ -b #X+)X+ + ] } , (14) 


where The vacuum pressure has been 

subtracted and at finite p the function <I> is generally dif¬ 
ferent from its complex conjugate d>. In the small den¬ 
sity limit constituent quark mass and expectation value 
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of traced Polyakov loop are ^ independent so the second 
derivative of the quark pressure simplifies giving after 
noticing that it can be written as derivative with respect 
to the quasiparticle energy Ep 

PE P 

^liT) = -g^^{T,^^ = 0) 


and the relation EpdEp = pdp has been used. Integrating 
by parts, 

ml{T) = r dppV^T^^^^ 

TT Jq dp 


SNftts 


pyjp^ + 711?/^ {Ep) 


r d^p d 

+ $(X -b 2A2) 


j (27r)3 dEp 

_l-b3$(A-bA2)-bX3_ 

, - J dpU(Ep)+ 


(15) 


where we have used X = and the fact that $ = $ 

for fi = 0. The quantity under the integral can be identi¬ 
fied with the energy derivative of our modified distribu¬ 
tion function bringing us to the following formula 


,(20) 


one arrives at the result for the Debye screening mass for 
Nf flavours of fermions 


I'D 


(T) = -c 


8N. 


f 


2 ^ , 

dpp ^(Ep) 


( 16 ) 


which after integration by parts (see section III) gives for 
Nf = 2 the result (11). 


,nl{T) = r dp [pi + E‘] . (21) 

^ Jo -^p 

For the non-abelian case the induced density reads 

Pi„d(a;) = ‘^9 J [fXiP.x) - ftip^x)] Tr[t^t“] , 

( 22 ) 


where Tr[f*’t“] = ^6°“^ and we assume 


III. KINETIC THEORY APPROACH 

The usual Debye screening mass in QED or perturba¬ 
tive QCD can be derived within a kinetic theory approach 
PSI [55] . In this section we will modify the standard ki¬ 
netic theory so that we consistently reproduce our previ¬ 
ous Debye mass derivation and generalize it in order to 
make contact with perturbative calculations for high tem¬ 
peratures. The kinetic theory approach has been widely 
used in the context of perturbative QCD Ha 130] pro¬ 
viding a physical picture behind the hard thermal loop 
approximation and some insights into transport prop¬ 
erties and collective modes of the quark gluon plasma. 
The appropriate change with respect to the textbook re¬ 
sult is that the usual Fermi Dirac distribution function 
is replaced by the Polyakov loop modified distribution 
function (12|. Then the charge density induced by the 
electrostatic potential Ao{x) = V{x) can be written 

Pi„d((i;) = 25 y [f’S-iEp + gV{x)) - U{Ep - gV{x))] 
2 / d^P dU 


f±{x,p) = ±g^{Ep)V‘^{x) . 


(23) 


Doing the same steps as before will give the Debye mass 
{Nf = 2) 




{T)=^ rdp[p‘+Ei] 

^ Jo 


MEp) 

En 


(24) 


which now reproduces the Debye mass of perturbative 
QCD for high temperatures. 


IV. COMPARISON WITH LATTICE QCD 

Within IQCD, the temperature dependent screening 
masses have been defined from the exponential fall-off of 
the colour singlet free energies [3] and the results could 
be represented in the rescaled leading order perturbative 
result 






(27r)3 dEp 

where the factor 2 is due to the fermion spin. The 
Maxwell equations give 

-V‘^V{x) = pi„d(a;) = ^ [ dpp^^{Ep)V{x) 

Jo dEp 


(25) 


= -ml{T)V{x) , 


(18) 


where the factor A{T) was introduced to account for non- 
perturbative effects. The analysis performed in m has 
shown that A{T) « 1.66 for Nf = 2-1-1 and T > 1.2 Tc- 
Obviously, A{T) —>■ 1 for high temperatures in agreement 
with perturbation theory. In order to compare our model 
with IQCD data we have to adopt equation (24), because 


where 




(T) = - 


—L 

EdppVp^^Al 09) 

TT Jo dp 


it makes contact with perturbation theory for high tem¬ 
peratures. The running coupling constant is modelled in 
two forms. The first one being the two-loop result [33] 


a{T) = 


1 


47r 


2/3oln(f ) + (f^)ln(21n(f)) 


(26) 
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where fj, = tt is the upper bound for the perturbative 
coupling. For A we use a value in the MS scheme of 
^Ws ~ MeV and 


/3o 


1 

IGtt^ 



(27) 


/3i 



38iV/ ^ 


(28) 


The other possibility is to use a running coupling con¬ 
stant obtained by solving the one-loop renormalization 
group equation with pole subtraction [55] 


as{T) 


47r 1 

Jo 21nf + A2-(^T)2 


(29) 


where /3o = 11 — 2/3Nf. 

In this way we can identify our model predictions for 
the nonperturbative effects of confinement and chiral 
symmetry restoration which are expressed as a tempera¬ 
ture dependent factor A{T) which reads (for Nf = 2) 

‘A + 

We see that it is entirely controlled by the quark distri¬ 
bution function and by the temperature dependent quark 
mass which mimic confinement as well as chiral symme¬ 
try breaking aspects of strong interaction dynamics. 

For a comparison to lattice QCD data we have chosen 
the Nt = 6 data for 2-|-1-flavors from [T5|. There two def¬ 
initions of the running coupling were used, as from a fit 
using a Debye screened Coloumb Ansatz at large separa¬ 
tions and amax obtained at the maximum of the effective 
r and T-dependent running coupling. Those results are 
based on an analysis of gauge field configurations gener¬ 
ated by the RBC-Bielefeld collaboration in (2-|-l)-flavor 
QCD for the calculation of the QCD equation of state [55] 
where the pion mass is about 220 MeV and the strange 
quark mass is adjusted to its physical value. 

From the lower panel of Fig. we see that the gen¬ 
eral trend of the lattice data is reproduced. The fact 
that our result for the Debye mass stays below the lat¬ 
tice data is due to the lack of dynamical gluons which 
in full QCD also bring a contribution to the screening 
effects. Also the smoother behaviour around is due 
to the used form of the running coupling which in prin¬ 
ciple is not applicable in the transition region. Below Tc 
the suppression of colour charges (which is meant as a 
rough form of confinement) drives mjo{T) fastly to zero 
and overcompensates the increasing interaction strength 
which taken alone would tend to increase the screening 
mass. Note that the quenched results for tod in [T2| and 
m lie below the lattice data shown in Fig. due to the 
missing degrees of freedom dynamical quarks. 




FIG. 2: Upper panel: Two loop running coupling an the an¬ 
alytical running coupling |34[ I35j . Lower panel: Debye mass 
with coupling to the Polyakov loop and running coupling con¬ 
stant a{T) compared with the lattice data [3]. 


V. CONCLUSIONS 

In this short note we have presented a calculation of 
the quark contribution to the Debye screening mass us¬ 
ing the PNJL model which captures chiral dynamics and 
in a very rough way some aspects of confinement. We 
have compared our results to the Debye masses extracted 
from lattice QCD simulations of the static heavy-quark 
potential and obtain overall agreement for the shape of 
the temperature dependence. Naturally, our results for 
the Debye mass stay below the lattice results since in our 
model the gluon contribution is neglected. However, as 
the observed gluon contribution on the lattice is of the 
same shape m we are convinced that our model cap¬ 
tures the essentials of the influence of chiral dynamics 
and confinement on the screening of the heavy-quark po¬ 
tential. A further improvement of the calculation would 
be to include dynamical gluons into the system, to im- 
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prove the modelling of quark (and gluon) confinement 
and to elaborate on the behaviour of the running cou¬ 
pling constant. The latter should also be compared with 
the lattice QCD result for this quantity as it is obtained 
simultaneously with that for the Debye mass. At this 
point it is interesting to go back to the different inter¬ 
pretations of the same lattice data. Here we would like 
to mention the one by Ref. m where the ansatz of a 
screened Corrnell type potential was adopted with two 
different Debye masses, one for the linear confining part 
(to£)) and one for the Coulombic part The per¬ 

formed fit gave a drastically different behaviour of the 
two screening masses. The temperature dependence ob¬ 
tained for fhu appears similar to that of the Debye mass 
in the present approach, calculated for a T- indepen¬ 
dent coupling constant (see Fig.[^. The physical reason 
for such a distinction could be that the stringy and the 
Coulombic parts of the potential act on different length 
scales so that the screening of them involves different dy¬ 
namics. The linear part should be dominant for larger 
distances thus involving stronger interactions and more 
correlations in the screening mechanism. Thus one could 
expect that rfiD > rno for all temperatures which is the 
finding of HU. This is, however, only a qualitative ar¬ 
gument. Our calculation, since it is at one-loop order 


with a QED like interaction should apply to the screened 
Coulomb potential part from which the lattice QCD re¬ 
sult has been extracted. As has been demonstrated here, 
this comparison provides a reasonable interpretation of 
the temperature dependence of the Debye mass. 
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Appendix A: Polarization loop, temporal component 


The calculation we have performed follows in all steps the standard QED evaluation of the polarization loop 
with the only difference that the usual Matsubara summation is now equipped with a trace over the colour indices 
|23| . To start with let us define the 44 component of the polarization tensor 

1144 = ^ j ^^^Tr| [ 74(771 + 74W„ - 7^74(m + 74a;„ - 74Wi - 7 (p-^)] , (Al) 

where Tr stands for the trace in Dirac and color spaces, a;„ = {2n + 1 )ttT — A 4 , with the temporal gluon held A 4 , 
uji = 2'kIT. Let us dehne 




1 


_ ‘2iEp 'lUJ'fi ~\~ sEp 


(A2) 


u;„ + 

where Ep = y/p'^ + m?. The hrst step is to calculate the Dirac trace using {7,^,7^} = — 25 py 

Tr [(m + 740;^ + jp)'yl{m + 740;^ - 740;; - j{p- ^)] = -4 [tti^ - a;„(w„ - w/) + p{p- g)] , (A 3 ) 


A/y = w„(a;„ - uji) - p{p- (p) - - cuntoi -p^ +pq- 


Tr [(tti + 74 W„ + pp)jl{m + 74 W„ - 74 a;/ - 7 (p - q))] = 4A/V • 


(A4) 

(A5) 


q'^T 

n44(*wj;q) = —T 




2A(iuj„,p) + {4pq - AEI - q^ - ujf)A(iuJn,p)A(iuJn - uJhP- q) > ■ (A 6 ) 


Now we want to decompose (A34) so that we have A{iu!mP) + A{iujn — iuji,p— q) and for that we introduce 

M =p'^ + 2nP + + (w„ - uJiY + [p- qf' = 2.{p^ + - UnUJi - pq) + wf + q^ , (A7) 


2Mv = Ai + Apq — Ap^ — ArrP — q^ — ojf . 
Further, we evaluate three coloured Matsubara sums 


J 3 = T^Trc [A( 7 a;„,^A( 7 w„ - iuji,p- q)] 


(A 8 ) 


(A9) 


Ji=T '^TicA{ioJn,p) ■ 


(AlO) 


J 2 = r^Trc{[p(p- ^]A( 7 a;„,^A(zw„ - iuji,p- q)} , 


(All) 


following [23] by going to the Polyakov gauge (where Polyakov loop variable is diagonal) and explicitly calculating 
the trace. The outcome reads 


^3= E 


55 


^^ ^EpEp—q tUJ “ 1 “ sEp 5 Ep—q 


[MsEp)-Ms'Ep_g)] , 


(A12) 


J 3 = 


AE E 

•-tj^pj^p—q 


[U{Ep) - U{Ep.,)] 


1 


1 


[I- U{Ep)- U{Ep_,)] , 


KjJ ~\~ Ep Ep — q ZUJ Ep Ep — q 

1 1 


+ 


iuj + Ep + Ep—q 


iuj Ep Ep—q 


(A13) 
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J2 = 


55 


p(p- q) 


^ , AE„Er,-n iu! + sEj, — s' E„-^ 


[U{sE^) - h{s'E^_^)] 


s,s'=± 


(A14) 


J 2 = f^ j IMEp) - MEp_g)] 


AE E 

^±^p±^p_q 


[1 - MEp) - MEp_,)] 


iu) + Ep — Ep_q iu) — Ep + Ep_q 


1 


+ 


%(jJ -\~ Ep Ep — q 


icj Ep Ep—q 


(A15) 


1 


*^4 = E = ^[1 - ^MEp)], 

^Ep ZEp 


(A16) 


which is in agreement with (6.27) and (6.37) from LeBellac, with the only difference that Fermi-Dirac distribution 
functions have been replaced with the modified ones [53] 


/$(A)=r^Tr, 


1 


iuJn — E 


= 3 


$(l + 2e-^4J)g-/3B_^g-3/3B 


1 + 3$(1 + e~P^)e~P^ + e“3/3£^ 


(A17) 


To obtain the last equation we use the fact that in this specific gauge the Polyakov loop variable is diagonal and that 
after a Matsubara summation we get 


U{E) = Y.- 




—' 1 + B dE 

=1 i=i 

where and A is to be understood as a temporal component of the gauge field, 

colour trace is now trivial and results in 


^ln(l + L^je '^^) = In [(1 + Liie ^^)(1 + A22e ^^)(1 + L33e ^^)] 
i=i 


(A18) 

The evaluation of the 

(A19) 


Using Lii + L 22 + T 33 = TrL = TrL^ = 3$, T 11 L 22 T 33 = det L = det = 1 we get 

h{E) = + 3<J>(1 + , 


which is the modihed distribution function (12). We also show that 

fq,{x + iuji) = f^{x), /$(-A) = 1 -/$(F;) . 


(A20) 


(A21) 


Now to get rid of the f(Ep-q) terms whenever we meet them we change the variables according to [31] p —)■ —p + q 
so that Ep_q — 7 > Ep, Ep Ep_q and p{p— q) does not change. We also make the Wick rotation ioj ^ oj 


J 3 


MEp) 

E„ 


+ 


(u; + Ep)^ - E^_q {uj-Ep)^-E^_q\ 


(A22) 


J 2 = p{p- ^ 


UjEp) 

Er, 


{u: + EpY-El_q {uj-EpY-El_q 


= p{p- ^Js 


(A23) 


E^ - El_q = 2pq -q^ pq= pqX 



dpp^dX 


(A24) 


(w ± EpY - Ep_q = Uj'^ - q^ ± 2ujEp + 2pqX , 


(A25) 
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2.pqX —y A . 


(A26) 


We are now left with two angular integrals 


f 2 pq 


l-2pq 


dX 


= In 


— ( 7 ^ ± 2uiEp + 2pq 


(^2 _ q ,2 2ujEp + X — ( 7 ^ ± 2ujEp — 2pq ’ 


(A27) 


f 2 pq 
- 2 pq ‘ 


XdX 


l- 2 pq a ;2 - ± 2 LoEp + X 


A r 2 2 _L o 11 + 2pq 

= 4pq - w -q ± 2a;£'p In ^ ^-l o f ^ 

— q^ ± 2u}Ep — 2pq 


(A28) 


R± (w) = — q^ — 2ujEp ± 2p(7 . 


(A29) 


Using the definition Ref{oj) = ^[/(w) + /(—w)] we obtain for the longitudinal component of Q the result 


F{u;,q) = -g^^^^^Re 


roc> ^2 


p'^dp2f,f,{Ep) j ^ ^ 4:EpOJ + q-^ - uj^ - AE^ 


En 


1 + 


Apq 


R_{uj) ■ 


(A30) 


To get the electric screening mass (Debye mass) we have to compute U(0, g —>■ 0) = mj^, i.e., 


F{ 0 ,q) = 9^ [ 

Jo 


°° p'^dp UiEp) \ q^-AEl {2p-qf 


- ^ ; 2 + 4_ii In — _ — \ 

^2 1 ^ Apq ( 2 p + g) 2 [’ 


(A31) 


q^o q ( 2 p + gj"= p 


(A32) 


I 7 ('n V n\ _ ™2 _ 2 f dp 2f^{Ep) ^ 2 , 17,21 
F{0,q ^ 0) — rrijj — g / — - [p + Epj . 

Jo ^ 


(A33) 


This result has a structure exactly as the QED Debye mass with the only difference that the Fermi-Dirac distribution 
function is replaced with the Polyakov loop suppressed distribution. After including the factors Nf = 2 and = g’^/Air 
we get 


16a, 


rrijj = 


dp 


f^{Ep) |.^2 I ^21 


E„ 


{P^ + E^p} , 


(A34) 


which is the formula 0 from the text. 

















